Tunneling time scale of under-the-barrier forerunners 
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Time-dependent analytical solutions to Schrodinger's equation with quantum shutter initial con- 
ditions are used to investigate the issue of the tunneling time of forerunners in rectangular potential 
barriers. By using a time-frequency analysis, we find the existence of a regime characterized by the 
opacity of the barrier, where the maximum peak of a forerunner measured at the barrier transmission 
edge x = L corresponds to a genuine tunneling process. The corresponding time scale represents 
the tunneling time of the forerunner through the classically forbidden region. 

PACS numbers: 03.65Bz,0365.Ca. 



Tunneling, that describes the possibility that a parti- 
cle traverses through a classically forbidden region, is one 
of the paradigms of quantum mechanics. In the energy 
domain one solves Schrodinger's equation at a fixed en- 
ergy E to obtain the probability of transmission through 
a barrier region, a subject discussed in every quantum 
mechanics textbook. In the time domain, however, there 
are still aspects open to scrutiny. A problem that has re- 
mained controversial over the years is the tunneling time 
problem, that may be stated by the question : How long 
does it take to a particle to traverse a classically forbidden 
region? Different authors have proposed and defended 
different views in answering the above question Q. Some 
authors have argued that it might not be a unique defini- 
tion of the tunneling time since different arrangements for 
the tunneling process may lead to different relevant time 
scales. An interesting idea that arises from the above 
considerations is that the relevance of the different tun- 
neling time definitions may depend on the context where 
they may become physically significant quantities^. 

In recent work we have investigated the effect of the 
transient solutions to the time-dependent Schrodinger's 
equation for cutoff wave initial conditions (quantum shut- 
ter) on the tunneling process^, ||, ||. In particular for 
tunneling through a barrier we found that just across 
the tunneling barrier, the probability density as a func- 
tion of time may exhibit a sharp peak maximum that 
we called time domain resonance. The peak value t max 
of that forerunner represents the largest probability of 
finding the particle at the barrier width L. We found 
also that the behavior of t max as a function of the bar- 
rier width L exhibits for small values of L a basin re- 
gion, followed by a region where t max grows linearly as 
L increases |I| . More recently, in collaboration with Del- 
gado and MugaQ, we have investigated the time scale 
for forerunners preceeding the main tunneling signal of a 
wave created by a source with a sharp onset j(| frj and by a 
quantum shutter B for systems with very large or infinite 
barrier width. For the particular case of opaque finite 



barriers we found a basin regime occurring at fixed po- 
sitions < x <C Kq 1 along the internal region, where fore- 
runners arc dominated by under-the-barrier frequency 
components. That is, in opaque barrier systems the tun- 
neling forerunners are observed at distances of the or- 
der or smaller than the penetration length Kq , where 
«o = [2m(V — E)] 1 / 2 /h, with E and V corresponding to 
the incidence energy and the barrier height, respectively. 
However, at distances k^ 1 < x < L the non-tunneling 
components eventually dominate the time evolution pro- 
cess, and hence it is not possible to speak of a genuine 
tunneling time scale of forerunners at x — L. We be- 
lieve that in order to gain more insight on the properties 
of transient tunneling structures and their corresponding 
time scales, the existence of tunneling forerunners at the 
barrier edge x = L needs to be investigated. 

The aim of this work is to show the existence of a 
regime along the basin region, where the maximum t max 
of the time domain resonance measured at the barrier 
width L corresponds to under-the-barrier tunneling and 
consequently provides a genuine time scale for tunneling. 

Our approach to the tunneling time problem is based 
on a model that deals with an explicit solution Q of time- 
dependent Schrodinger equation for an arbitrary poten- 
tial V(x) (0 < x < L) that vanishes outside the internal 
region. We consider the problem of the time evolution of 



a cutoff plane wave ^(x, k; t = 0) = 8(— x)(e 
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following the instantaneous opening at t = of a quan- 
tum shutter at x — 0. Along the tunneling region the 
solution reads, 

¥(x, k, t) = <!> k (x)M(y k ) - $_ k (x)M(y_ k ) 

OO 

- ]T $ n (x)M(y kn ), (0<x<L) (1) 
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where the ( I > ±fc(a;)'s refer to the stationary solutions of 
the problem and < &„(a;) = 2iku n (0)u n (x) / (k 2 — k 2 ). Sim- 
ilarly, the solution ty e (x, k; t) for the external or transmit- 
ted region (x > L), is given by B, 
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f e (x, k; t) = T k M(y k ) - T_ k M(y_ k ) 

OO 

-i T nM(y k „), (x>L) (2) 



where the T± k , s refer to the transmission amplitudes, 
and the factor T n = 2iku n (0)u n (L) exp(— ik n L)/(k 2 — 
k 2 ). In Eqs. ([[]) and (S) the coefhcients $ n and T n are 
given in terms of the resonant cigcnfunctions {u n (x)} 
with complex energy eigenvalues E n — ti 2 k 2 l /2m, with 
(a n ,b n > 0). The resonant sums in Eqs. 
(|I|) and(Q) run over the full set of complex poles {k n }. 
The M's are defined as H, 
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where w is the complex error function [g] defined as 
w(z) = cxp(—z 2 )erfc(—iz), with arguments y q (x,t) = 
e~' l7T / 4 (m/2ht) 1 / 2 {x — hqt/m], where q = ±k, k± n . 

We shall explore the issue of the tunneling time scale 
in typical one dimensional potential barriers |j| of height 
V and thickness L, defined along the interval < x < L. 
Let us analyze the behavior of the maximum of the time 
domain resonance, t max , as a function of the potential 
barrier width, L. Let us recall that this transient struc- 
ture corresponds to the first maximum of the probability 
density I* 6 ] 2 , measured at the fixed position, x = L. 
This can be appreciated in Fig. [l] where we plot the 
L— dependence of t max (full dot), corresponding to po- 
tential barrier systems with parameters: V = 0.3 eV, 
incidence energy E = h 2 k 2 /2m = 0.001 eV, and effective 
mass for the electron m = 0.067m e . Here we can clearly 
observe a basin corresponding to a range of values of the 
barrier width. We can also appreciate that if L is fur- 
ther increased, t p grows linearly with L. Such a linear 
regime occurs at large barrier widths. We have argued 
j| that in this case, the tunneling process is inhibited 
and that the particle goes mainly over the barrier. We 
have recently suggested || that for small values of the 
barrier width L, the basin exhibited by t max , is a result 
of a subtle interplay between tunneling and top-barrier 
resonant processes. In what follows we shall investigate 
under what conditions the time scales associated to the 
basin, are in fact related to a tunneling process. We begin 
our analysis by choosing the case L = 4.0 nm depicted 
in Fig. [l], which is located around the minimum of the 
basin. In Fig. |^ we plot the normalized probability den- 
sity |^ e | 2 (solid line), normalized to the transmission co- 
efficient \T k \ 2 , as a function of time t for a fixed value of 
position x = L. We can clearly appreciate a time domain 
resonance peaked at the value t max — 5.17 fs. In order 
to determine if this structure is related to a genuine tun- 
neling process we analyze the frequency content of |4' e | 2 . 



FIG. 1: Maximum of the time domain resonance t max (full 
dot) as a function of the barrier width L, for an incidence 
energy E — .001 eV. In this case the barrier height is V = 0.3 
eV. See text. 
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FIG. 2: Time evolution of |* e | 2 (solid line), at x = L = 4.0 
nm. A full square indicates the position of the maximum of 
the time-domain resonance, at t ma x = 5.17 fs. 



We compute the local average frequency uj av j|, 10 
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and the instantaneous bandwidth o-]l(|, 
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where s = i, e according to consider, respectively, the 
internal or external solutions. 

In Fig. |^ we plot the relative average local frequency 
(relative frequency for short) uj av /ujy, where u>v = V/h 
is the cut-off frequency, along the relevant time interval, 



3 




Time (fs) Distance (x/L) 



FIG. 3: Relative average local frequency cu a v/^v (solid 
line) for the case depicted in Fig. ^. The cutoff-frequency 
lov jiov = 1 (dashed line) is included for comparison. In the 
inset we plot the instantaneous bandwidth a of the spectro- 
gram depicted in the main graph. Notice that the frequency 
deviations from t max are exactly zero, i.e., a(t max ) = 0. In 
all cases a full square indicates the position of tmax- 

discussed in Fig. ||. We can appreciate that in the vicin- 
ity of the maximum of the time-domain resonance, t max , 
the probability density is composed entirely by under- 
the-barrier frequency components i.e u> av /u)v < 1. This 
also occurs at the exact value t max , also indicated in the 
figure by a solid square. In the inset of Fig. || we plot the 
instantaneous bandwidth a of the spectrogram. Notice 
the absence of a frequency dispersion around the maxi- 
mum t max , i.e., <j(t max ) = 0. The above result indicates 
that in this case the peak of the time domain resonance, 
and the values close to it, refer to a genuine tunneling 
event. In Fig. [| we plot the relative frequency uj av /ujv 
associated to different values of the maximum t max of the 
time-domain resonance, measured at different positions 
along both the internal and external regions of a poten- 
tial barrier with parameters: V — 0.3 eV, L = 4.13 nm, 
and a = 3.0. In this case we choose the following values 
of the incidence energy: E = 0.001 eV (solid dot), and 
E = 0.01 eV (hollow dot). In the inset of Fig. || we show, 
for the particular case of E = 0.01 eV, the values of t max 
(solid square) at the different values of position consid- 
ered in the main graph. As can be clearly appreciated in 
that figure, the tunneling process along the whole inter- 
nal region is governed by under-barrier-frequency com- 
ponents, i.e., uj av /ujy < 1. The above results should be 
contrasted with those obtained for very large or infinite 
barrier widths recently studied in Ref. where under- 
the-barrier tunneling components arise only at distances 
smaller than a characteristic length. Here we have found 
a particular combination of potential parameters char- 
acterizing a genuine tunneling process. The existence of 
under-the-barrier tunneling forerunners is not only re- 



FIG. 4: Relative frequency u a v/^>v of the maximum of the 
time-domain resonance, as a function of position, measured in 
units of the barrier width L. The parameters are given in the 
text. Two incidence energies are considered: E = 0.001 eV 
(solid dot), and E — 0.01 eV (hollow dot). In both cases, the 
relative frequency u a v/^v along the internal region is below 
the cutoff-frequency u>v l^v = 1 (dashed line) . The behavior 
of t max as a function of position is illustrated in the inset for 
the case with E — 0.01 eV. The position of the barrier edge, 
x — L, is indicated by a dotted line in both figures. 
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FIG. 5: Relative frequency u a v/^>v measured at the barrier 
a function of the opacity a. Here we considered 
a barrier height V = 0.3 eV, and the parameters: u = 300 
(solid dot), u — 10 (solid triangle), and u = 5 (solid square). 
Note that for values of the opacity smaller than a ~ 3.3, the 
relative frequencies for all values of u are below the cutoff- 
frequency ujv/lov = l(dashed line). See text. 



stricted to values of the position along the internal re- 
gion of the potential. In fact, we can see in Fig.||, that 
for distances up to x ~ 2L along the external region, we 
can still observe frequency components below the cutoff- 
frequency, ojv- As x/L increases further, Lj av /u)v > 1. 
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This behavior follows because the corresponding values of 
tmax increase linearly, indicating the prevalence of non- 
tunneling components (see the inset in Fig. ^). 

We have found that the existence of tunneling forerun- 
ners at the barrier edge x — L may be described more 
generally by referring to the opacity a of the system, 
defined as, 



[2mV] 



1/2 



-L. 



(6) 



and by the dimensionless parameter u, the ratio between 
the potential barrier height and the incidence energy, 



V 
E' 



(7) 



There is a regime, that we shall refer to as tunneling 
regime, characterized by the opacity a and the parameter 
u, where the relative frequencies associated to the time 
domain resonance are below the cutoff- frequency wy. In 
order to show this, in Fig. || we plot the relative fre- 
quency uj av /ujv as function of the opacity a, for three 
different values of the parameter u, u = 5 (solid dot), 
u = 10 (solid triangle), and u = 300 (solid square). Al- 
though in this case we have chosen a value of V = 0.3 
eV, we have checked that the relative frequency u} av /ujy, 
exhibits a striking regularity when plotted as a function 
of a. For a given value of a all the systems with the 
same parameter u, share the same curve. That is, all 
systems characterized by the same parameters a and u 
yield the same relative frequency, Lu av /ujy. In view of 
this observed regularity, we can completely characterize 
the regime associated with under-the-barrier frequency 
components. Although in Fig. |^ we have considered val- 
ues of the parameter u such that 5 < u < 300, the case 
corresponding to very large values of u (it — > oo), almost 
overlaps with the case u = 300. Thus, for very large val- 
ues of u there is an upper-bound for the opacity a u ~ 3.3. 
This result, and the fact that the lower-bound for the 
opacity [1TJ is given by a c — 2.065, recently reported 
in Ref. [ 5[ , allows to characterize an opacity "window" 
given by a c < a < a u where the relative frequencies are 
always below the cutoff- frequency ujy /toy — 1, irrespec- 
tive of the value of the parameter u, namely of the value 
of the incidence energy. Note that the values of a within 
the opacity "window" may be obtained using typical pa- 
rameters of semiconductor heterostructures^J. From the 
figure one can appreciate the existence of another regime 
of a's where the existence of under-the-barrier processes 
associated to the time-domain resonance may depend on 
the value of the parameter it. 

It is of interest to remark that the time scale given by 
the peak maximum t max of the forerunner at the barrier 
edge x = L, i.e., the time domain resonance, differs from 
both the semi-classical Biittiker-Landauer and Biittiker 
traversal times, which exhibit a linear dependence with 



L, and also differs from the phase-time^, ||]. In Ref. j|] 
we have discussed an unexpected relationship between 
the delay time and the existence of time domain reso- 
nances. 

To conclude we remark that the analytical solution to 
the time-dependent Schrodinger's equation with quan- 
tum shutter initial conditions applies in general to arbi- 
trary potentials provided they vanish beyond a distance 
and can also be extended to deal with cutoff pulses as 
discussed in Ref. || . Our results refer to novel transient 
effects in time-dependent tunneling. To test our results 
experimentally would require to consider the detection of 
tunneling particles in time domain at distances close to 
the interaction region. 
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